The Gouy phase shift in nonlinear interactions of waves 
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We theoretically analyze the influence of the Gouy phase shift on the nonlinear interaction between 
waves of different frequencies. We focus on \^ interaction of optical fields, e.g. through birefringent 
crystals, and show that focussing, stronger than suggested by the Boyd-Kleinman factor, can further 
improve nonlinear processes. An increased value of 3.32 for the optimal focussing parameter for a 
single pass process is found. The new value builds on the compensation of the Gouy phase shift by 
a spatially varying, instead constant, wave vector phase mismatch. We analyze the single-ended, 
singly resonant standing wave nonlinear cavity and show that in this case the Gouy phase shift leads 
to an additional phase during backreflection. Our numerical simulations may explain ill-understood 
experimental observations in such devices. 
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Nonlinear interactions of waves, in particular those of 
optical fields, have opened new research areas and have 
found various applications. In general, waves of different 
frequencies are coupled via nonlinear media, like bire- 
fringent crystals. Examples are the production of higher 
harmonics of laser radiation [l[ , the generation of tunable 
frequencies through optical parametric oscillation Q and 
the generation of nonclassical light 0, |j| for high preci- 
sion metrology 0, Q , fundamental tests of quantum me- 
chanics and quantum information [3] • The efficiency of a 
nonlinear process depends on parameters of the nonlinear 
medium, and generally increases with higher intensities 
of the fields involved and with better phase matching of 
their wave fronts. To achieve strong nonlinear interac- 
tions, pulsed laser radiation, strong focussing and, espe- 
cially for continuous wave radiation, intensity build-up 
in resonators are used. In plane wave theory perfect 
phase matching is achieved if the wave fronts of inter- 
acting fields propagate with the same velocity. For fo- 
cussed laser beams, however, this is not true because of 
the well-known Gouy phase shift. This phase shift oc- 
curs due to the spatial confinement of a focussed wave 
and generally depends on the spatial mode as well as the 
frequency of the wave Q . The influence of focussing into 
a nonlinear medium has been investigated by Boyd and 
Kleinman in great detail [jj]. They discovered that the 
efficiency of the nonlinear process does not monotonically 
increase with decreasing focal size. They especially con- 
sidered the lowest order nonlinearity that enables second 
harmonic generation (SHG) and optical parametric am- 
plification (OPA) and is described by the susceptibility 
X^ 2 \ and numerically found an optimum factor between 
the length of the nonlinear crystal and the Rayleigh range 
of the focussed Gaussian beam for a single pass through 
the crystal. 

In this Letter we show that focussing stronger than 
suggested by the Boyd Kleinman factor can further im- 
prove nonlinear processes. We show that this effect can 



be understood by considering the Gouy phase shift be- 
tween the interacting waves. We also show that the Gouy 
phase shift results in a non-trivial phase mismatch prob- 
lem in standing wave cavities. 

Boyd and Kleinman have found that the maximum 
nonlinear coupling between two Gaussian beams of fun- 
damental (subscript 1) and second harmonic waves (sub- 
script 2) is achieved for a positive wave vector phase mis- 
match Afc = 2|fci| — 1&2 >0 which increases with decreas- 
ing waist size of the beam. For a single pass through a 
nonlinear crystal of length L they numerically found the 
optimal focussing parameter given by the relation 
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where zr = Trw^n/X is the Rayleigh range of the beams 
inside the crystal, and wq, n and A are the beam's waist 
size, refractive index and wavelength, respectively We 
first show that the Boyd-Kleinman factor according to 
Eq. (JTJ) is a consequence of maximizing the intensity of 
the mean pump field inside the nonlinear medium under 
the constraint of the Gouy phase shift. In a medium 
the nonlinear interaction is described by the following set 
of differential equations 
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where i?o,i, ^0,2 are the electrical fields of the fundamen- 
tal and the harmonic mode in the focal center at position 
Zq, and Afc = Air/XAn is the phase mismatch between 
the two interacting modes. Here we use the following 
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abbreviations 



A(fr(z) = — arctan 



Z - Zq 

zr 



w(z) = w \ 1 + 
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where w(z) corresponds to the beam width at the posi- 
tion z. In plane wave theory one finds g{z) = cxp(iAkz) 
and A(fr = 0, and equal indices of refraction for the 
two interacting modes provide the maximum nonlincar- 
ity. When focussing the beam into a nonlinear material 
however, there is a non zero phase difference A(fr. Since 
the phase difference between a plane wave and a focussed 
Gaussian beam is given by the Gouy phase shift, A(fr 
should be the difference of two such phase shifts. When 
considering phase shifts between different oscillator fre- 
quencies, phases have to be frequency normalized. We 
therefore introduce the Gouy phase shift normalized to 
the optical frequency of mode i 
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where m and n describe the spatial Hermite-Gaussian 
modes (TEM„ m ). If we now normalize A(fr to the har- 
monic frequency we find 



A(j> 
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We point out, that for an optimized nonlinear interaction 
of Gaussian beams the Raylcigh ranges are identical for 
all modes involved. In the case of frequency conversion 
of a single pump field this is automatically realized by 
the nonlinear process. For the x -processes considered 
in Eqs. ©-(g]) (m = n = 0) we find 



A<t> = <j>G • 



(9) 



From this one can conclude that the Gouy phase shift 
leads to a nonperfect matching of the (nonplanar) phase 
fronts in nonlinear processes. To quantify this effect we 
define the effective nonlinearity k of the process. This 
quantity is proportional to the conversion efficiency in 
SHG as well as to the optical gain of OPA. For weak 
interaction, i. e. the pump field is not depleted by the 
nonlinear interaction, k is given by 
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where the integration is taken over the whole interaction 
length. For a single pass through a nonlinear medium of 
length L the effective nonlinearity is given by 
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FIG. 1: (Color online) Left: For weak focussing into the non- 
linear medium (£ = 0.18) the Gouy phase shift can be com- 
pensated by choosing AA: = 1/zr, and perfect phase match- 
ing can be realized over the full crystal length. Right: For 
stronger focussing (£ = 2.03) a constant Ak can not provide 
perfect phase matching, (a) Gouy phase shift A(fr, (b) com- 
pensating phase Akz, (c) overall phase (fro, where a constant 
value describes perfect phase matching. 



This quantity is maximized if the averaged field 
strength inside the crystal is maximized, i. e. if the focus 
is placed in the crystal center and if the condition 



Akz + Aip(z) = (fro = const. 



(12) 



is satisfied. In this case all partial waves are produced ex- 
actly in phase to each other, and perfect phase matching 
is realized. 

Curves (a) in Fig. Q] show the differential Gouy phase 
shifts Acfr(z) = 4>g{z) for weak and strong focussing, 
respectively. For weak focussing the gouy phase shift 
evolves linearly inside the medium, and one can compen- 
sate this phase mismatch by choosing Ak = 1/zr > 0, as 
found by Boyd and Kleinman For stronger focussing, 
however, it is not possible to achieve perfect compensa- 
tion from Afc that is constant over the crystal. Curves (b) 
show the compensating linear phase Akz that is due to 
the propagation inside the medium and curves (c) show 
the total phase cV The value of Ak was chosen to pro- 
vide the lowest variance of (fro over the whole interaction 
range. 

We now show that it is possible to realize perfect phase 
matching for an arbitrary focussing by applying the fol- 
lowing position dependent index of refraction 
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where the constant value of the phase (fro is set by the 
Gouy phase at the entrance surface of the nonlinear 
medium 
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Fig. [2] shows An sp (z) for a nonlinear crystal of length 
L for different focussing parameters £, with focus placed 
into the crystal's center. Curves (a) to (c) could experi- 
mentally be realized by applying an appropriate tempera- 
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FIG. 2: (Color online) Change of refractive index along the 
crystal to compensate for the Gouy phase shift, for three dif- 
ferent strengths of focussing. The temperature scale on the 
right corresponds to a MgO(7%) : LiNbOs crystal that has 
been used in @|. (a) £ = 0.55, (b) £ = 1.14, (c) £ = 3.32. 



ture gradient along z-dircction. Alternatively, an electri- 
cal field applied to the crystal could be used. The tem- 
perature values on the right vertical axis constitute an 
example for 7 % magnesium-oxidc-doped lithium niobate 
(MgOiLiNbOs). Data was deduced from measurements 
of the nonlinear efficiency of crystals used in type I OPA 
in @]. 

The effective nonlinearity for the single pass setup with 
perfect Gouy phase compensation can be written as 
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with £ = L/2zr. The numerical optimization of the 
above expression leads to 
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opt 



3.32. 
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Our result means that with optimum, position dependent 
phase matching, the optimal waist size is approximatly 
7.5 % smaller than suggested by the Boyd-Klcinman fac- 
tor, and according to this, the effective nonlinearity is 
further increased by 4.4%. 

We now analyze if the nonlinear interaction in 
standing-wave cavities can be similarily improved. 
Standing wave cavities, in particular in the form of a 
singly-resonant, single-ended cavity, i. e. with one mir- 
ror of almost perfect reflectivity, are frequently used in 
quantum and nonlinear optics Jg, [l(j ■ In such cavities 
waves that propagate in two different directions interfere 
with each other, and the differential phases introduced 
by the reflections at the cavity mirrors have to be con- 
sidered. Paschotta et al. have investigated the phase 
difference Aip introduced from back reflection and have 
suggested an appropriate design of the high reflectivity 



dielectric multi-layer coating to annihilate any additional 
phase shift. It can be shown that the effective nonlinear- 
ity for a doublepass of plane waves through a crystal 
depends on the Aip in the following way 
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For the calculation of the doublepass effective nonlinear- 
ity in the case of focussed Gaussian beams we model the 
system with a nonlinear medium of length 2L with a thin 
lense at position L that refocuses the beam. In this way 
we obtain two waists at positions zq and z' Q = 2L — zq of 
size wo indicating the way to the endmirror and the way 
back, respectively. Now we integrate over 2L and find 
the following expression for the effective nonlinearity for 
a double pass of the fundamental Hermite-Gauss mode 
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where w'{z) and 4>' G {z) belong to the focus at position 
z' . Aip is again the differential phase that may be in- 
troduced by the coating of the back reflecting mirror. 
We first consider the special case of weak focussing, i. e. 
\z — zoI/zr, <§C 1 Vz G [0, L], and simplify the above 
expression as follows 
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where Afc' := Afc — 1/zr and 

Aip' := Aip + 2(L - z )/zr 



(19) 



(20) 



This expression has the same form as the one for plane 
waves as given in Eq. (|17p . However, an additional phase 
shift appears. This phase shift is a result of spatial con- 
finement and the swapping in sign of the wave front's 
radius of curvature during reflection, and corresponds to 
minus twice the Gouy phase in the limit considered here. 
From the expression of Aip' in Eq. (f20f it follows that this 
additional phase jump vanishes if the waist is located ex- 
actly at the back reflecting surface. In this case we have 
plane wave fronts at the end mirror and therefore the 
system is similar to a single pass through a nonlinear 
medium of length 2L. 
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FIG. 3: (Color online) Effective nonlinearity Kd P (normal- 
ized) versus the differential phase introduced by the back re- 
flecting surface Aip. For Gaussian beams with waist position 
on the reflecting surface, Aip = generally provides the high- 
est Kdp, ( a ) (here with strong focussing £ = 2.84, zo = L); the 
same result is found for plane waves. Contrary, for £ = 2.84 
and zo — L/2 we find Aip = n (c). The result for weaker 
focussing is shown in (b) (£ = 0.775, zo = L/2). 



We now examine Kdp for reflected Gaussian beams of 
arbitrary focussing parameters £ and the waist located in 
the center of the medium. Eq. (|18|) yields 



^dp = 4 cos^l ■ 
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dz 
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The first term provides the optimal phase of the end mir- 
ror of Aip = —AkL. In turn Afc is again found by mini- 
mizing the variance of the term Akz + 4>g( z )- We obtain 
the following expression for the optimal differential phase 
Aip of the end mirror versus focussing parameter 

A^ = - A [(1 + e) ■ arctan(C) - £] . (22) 

Fig.[3]shows the effective nonlinearity versus Aip for three 
different standing wave cavity arrangements. In all cases 
the second harmonic wave is not resonant but simply 
back reflected. Curves (a) and (c) use the focussing pa- 
rameter £ = 2.84. This value optimizes the effective 
nonlinearity of the cavity if the refractive index of the 
medium does not depend on direction of propagation of 
waves. This is evident from Fig. [2] because the position 
dependent refractive indices are not symmetric with re- 
spect to the focal position and the back reflected wave 
would require different values. However, if one transfers 
the results from a single pass through the medium and 
realizes the required propagation direction dependent re- 
fractive index the optimum focussing parameter is again 
£ = 3.32. Curve (a) represents the case for focussing di- 
rectly onto the back reflecting surface. Curve (c) shows 



the effective nonlinearity for a waist position at the crys- 
tal's centre. In the latter case the best choice of the dif- 
ferential phase at the back reflecting surface is Aip « n. 
This is exactly the opposite of what one might expect 
from plane wave theory, where the optimum phase is 
Aip = 0, similar to curve (a). Trace (6) shows the ef- 
fective nonlinearity for the focussing parameter that was 
chosen by Paschotta et. al. [llj . In that paper a full 



quantitative comparison between experiment and theory 
of the nonlinearity in standing wave cavities was con- 
ducted. In their experiment the back reflecting mirror 
was designed to prevent a differential phase shift between 
the two interacting modes and a value of Aip = was 
chosen. However, their experimental data revealed the 
effective nonlinearity to be « 10 % smaller than expected 
from their calculations. From our calculation it follows 
that the optimum phase for their setup was Aip w 1.55 n 
and that the chosen value of Aip = decreased the effec- 
tive nonlinearity to about 90 % of the maximum value. 
Our considerations are therefore in excellent agreement 
with experimental results given in that paper and can 
solve the observed discrepancy 

In conclusion we have shown how for focussed waves 
the Gouy phase shift produces nonideal phase matching 
in case of Afc = 0. For a single pass through a nonlin- 
ear medium the optimum focussing parameter is found 
to be £ = 3.32. In this case a position dependent re- 
fractive index is required to further improve the effective 
nonlinearity by 4.4%. For a double pass and for cav- 
ities an optimum focussing parameter above £ = 2.84 
can only be achieved with a refractive index that also de- 
pends on propagation direction. We have also shown that 
the Gouy phase shift effects the optimum value for the 
phases introduced by cavity mirrors, with a significant 
effect on the achievable effective nonlinearity. Our theo- 
retical analysis shows exact agreement with experimen- 
tal data published elsewhere, and may lead to improved 
quantitative descriptions of nonlinear cavities. 
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